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Abstract. The time-evolution of the maximum and the width of exact analytic wave packet 
(WP) solutions of the time-dependent Schrodinger equation (SE) represents the particle and 
wave aspects, respectively, of the quantum system. The dynamics of the maximum, located 
at the mean value of position, is governed by the Newtonian equation of the corresponding 
classical problem. The width, which is directly proportional to the position uncertainty, 
obeys a complex nonlinear Riccati equation which can be transformed into a real nonlinear 
Ermakov equation. The coupled pair of these equations yields a dynamical invariant which 
plays a key role in our investigation. It can be expressed in terms of a complex variable that 
linearizes the Riccati equation. This variable also provides the time-dependent parameters 
that characterize the Green's function, or Feynman kernel, of the corresponding problem. 
From there, also the relation between the classical and quantum dynamics of the systems 
can be obtained. Furthermore, the close connection between the Ermakov invariant and the 
Wigner function will be shown. Factorization of the dynamical invariant allows for compari- 
son with creation/annihilation operators and supersymmetry where the partner potentials 
fulfil (real) Riccati equations. This provides the link to a nonlinear formulation of time- 
independent quantum mechanics in terms of an Ermakov equation for the amplitude of the 
stationary state wave functions combined with a conservation law. Comparison with SUSY 
and the time-dependent problems concludes our analysis. 

Key words: Riccati equation; Ermakov invariant; wave packet dynamics; nonlinear quantum 
mechanics 
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1 Introduction 

Participating in a conference on nonlinear mathematical physics taking place in Kiev, one can 
find more colleagues than usual who are familiar with the name Ermakov; some are even able to 
tell where one can find traces of him from the time he used to work at the university in Kiev. This 
period was, however, before quantum mechanics was established and Vasili Petrovich Ermakov 
(1845-1922) died already four years before Schrodinger published his celebrated equation pQ. So, 
what can be his relation with quantum theory? The second name in the title of this paper also 
seems to be surprising in this context because Jacopo Riccati (1676-1754) was born actually 250 
years prior to Schrodinger's publication! The surprising answer to our question is that Ermakov's 
and Riccati's contributions are not only closely connected with quantum theory in the form 
established by Schrodinger and Heisenberg, but also with later formulations due to Feynman 
and Wigner. Therefore, their work plays a central role in finding underlying connections amongst 
all these different formulations. 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http://www.emis.de/journals/SIGMA/symmetry2007.html| 
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One way of showing how Riccati and Ermakov enter the formalism of quantum theory is the 
study of cases where exact analytic Gaussian wave packet (WP) solutions of the time-dependent 
Schrodinger equation (SE) exist, in particular, the harmonic oscillator (HO) and the free mo- 
tion. Therefore, in the following, we will restrict the discussion of the time-dependent problems 
to the (one-dimensional) HO (where the frequency uj may also be time-dependent) and the free 
motion is obtained in the limit u> — > 0. In Section 2, it will be shown how a (complex) nonlinear 
Riccati equation that governs the dynamics of a typical quantum effect, namely, the position 
uncertainty, emerges from the WP solution of the time-dependent SE. This Riccati equation 
can be transformed into a (real) nonlinear Ermakov equation that, together with the classical 
Newtonian equation of the system, provides a corresponding dynamical invariant, the so-called 
Ermakov invariant, that will play an important role in our discussion. As shown in Section 3, 
the linearization of the Riccati equation can be accomplished by introducing new complex vari- 
ables. These variables are connected with the time-dependent parameters that characterize the 
time-dependent Green's function, also called Feynman kernel, of the corresponding problem. 
They also provide the key to understanding the relations between the dynamics of the classical 
and quantum aspects of the system; in particular, the importance of the initial position un- 
certainty will become obvious. Furthermore, in Section 4, the close relationship between the 
Ermakov invariant and the Wigner function will be shown. The link to the occurrence of Ric- 
cati and Ermakov equations in time-independent quantum mechanics is obtained in Section 5. 
Factorization of the Ermakov invariant allows for a comparison with creation and annihilation 
operators as well as their generalizations used in supersymmetry (SUSY). In the latter case, the 
supersymmetric partner potentials fulfil real Riccati equations. This shows similarities with 
a nonlinear formulation of time-independent quantum mechanics by Reinisch [2], where he ob- 
tains an Ermakov equation (with spatial derivatives instead of temporal ones) for the amplitude 
of the stationary state wave functions. This Ermakov equation is related with a kind of com- 
plexification of the Riccati equation appearing in SUSY. Finally, in Section 6, the main results 
will be summarized and possible generalizations mentioned. 



2 Dynamics of Gaussian wave packets and Ermakov invariant 

Starting point of our investigation is the time-dependent SE for the HO with frequency u>(t), 



8 , f h 2 d 2 



m 



that possesses Gaussian WP solutions of the type 



% P (x,t) = N(t) exp <^ i 



y{t)x 2 + ^Lx + K{t) 
n 



(2) 



where y(t) = ypt(t) + iyi(t), x = x — (x) = x — rj(t) (i.e., (x) = rj(t) is the classical trajectory 
calculated as mean value J f^ p x$wp^ = (%))> {p) — m V is the classical momentum and N(t) 
and K(t) are purely time-dependent terms that are not relevant for the following discussion. 

Inserting the WP ([2]) into equation ([1]) yields the equations of motion for rj(t) and y(t). The 
equation for the WP maximum, located at x = rj(t), is just the classical equation of motion 

r} + u J 2 (t)r ! = 0, (3) 

where overdots denote time derivatives. 

The equation of motion for the complex quantity y(t) that is connected with the WP width 
and, thus, the position uncertainty, is given by the complex Riccati-type equation 

—y+[—y) +u 2 (t) = o. 4 



Riccati and Ermakov Equations in Quantum Theory 



3 



Equation Q can be separated into real and imaginary parts, 

2h [2h \ [2h \ 

Im : — yi + 2 — m —y R = 0, 5 
m \m J \m J 

2h. [2h \ 2 [2h 



Re : — y R + — y R — I — y\ I + u\t) = 0. (6) 
m \m J \m J 

The real part y R (t) can be eliminated from equation ([6]) by solving equation ([5]) for y R and 
inserting the result into equation ([6]). 

It is useful to introduce a new variable a(t) that is connected with yi{t) via 

m a z (t) 



1/2 



2m{x 2 {t)) 

h 



where a(t) is directly proportional to the WP width, or position uncertainty, i.e., a = ^2^7) 

1/2 

with (x 2 ) = (x 2 ) — (x) 2 . Inserting this definition (|7|) into equation ([5]) shows that 
the real part of y(t) just describes the relative change in time of the WP width, 

2h _a_li(x 2 ) 

m yR -«-2l^T (8) 
Together with definition ([7]), this finally turns equation (0) into 

a + oj 2 (t)a = 4r, (9) 
a 6 

the so-called Ermakov equation. 

It has been shown by Ermakov in 1880 [3] that the system of differential equations (|3|) and ([9]), 
coupled via the possibly time-dependent frequency uj, leads to a dynamical invariant that has 
been rediscovered by several authors in the 20th century [1], 



4-3 



(rja - -qa) 2 + ( — 



const. (10) 



It is straightforward to show that J^/l = 0; a proof following Ermakov's method can be found 
in [5]. The Ermakov invariant not only depends on the classical variables rj(t) and r](t), but also 
on the quantum uncertainty connected with a(t) and a{t). Additional interesting insight into 
the relation between the variables rj and a can be obtained from a different treatment of the 
Riccati equation (j4]). 

3 Linearization of the complex Riccati equation, 
Feynman kernel and quantum-classical connection 

Introducing a new complex variable X(t), the complex variable in the Riccati equation (j4|) can 
be replaced by the logarithmic time-derivative of A, i.e., 

— y )= t, ( n ) 

m J A 

thus turning the nonlinear Riccati equation into the complex linear equation 
\ + uj 2 {t)\ = 0, 

looking exactly like equation ([3|) for r](t), which is not just by accident, as will be shown later. 
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The complex variable X(t) can be written in polar coordinates as well as in cartesian coordi- 
nates, i.e., 



A = ae 



up 



u + iz. 



The choice of the symbol a for the absolute value of A is also not coincidental, as will now be 
shown. 

3.1 Polar form of the complex linearization variable 

Writing relation ([lip in polar coordinates yields 



A a 
A a 



+ i<p 



2h 



-1J 



m 



(12) 



The real part already looks identical to the one given in equation (jHJ). So, the absolute value of A 
would be (up to a constant factor) identical with the square root of the position uncertainty (x 2 ), 
if 



^ = ~2 
or 



(13) 



is fulfilled. This can easily be verified by inserting and (f^yi), as given in equation (fT2|) . 

into the imaginary part of the Riccati equation, thus turning equation ((SJ) into 



tp + 2-<p = 0, 
a 

in agreement with equation ([13j) . From equation for the real part then, again, the Ermakov 
equation is obtained as 



a + oj (t)a = atp 



•2 



1 



a- 



3.2 Cartesian form and Feynman kernel 

After the physical meaning of the absolute value of A in polar coordinates and its relation to 
the phase angle via (j!3p have been clarified, the interpretation of the cartesian coordinates u 
and z needs to be ascertained. For this purpose, it can be utilized that the WP solution ([2]) at 
time t can also be obtained with the help of an initial WP at, e.g., t' = and a time-dependent 
Green's function, also called time-propagator or Feynman kernel, via 



*wp(x,t) = J dx'G(x,x',t,t' = 0)^ W P (a/,0) . 
For the considered Gaussian WP with initial distribution 



$wp (x',0) 



exp 



im 
2h I 



■a /2 i n P° I 

ipox + 2 — x 
m 



(14) 



(15) 



where (3q = 2m (^) = ^ an< ^ Po = (p){P = 0)> the Green's function can be written as 



im 



exp <( — 



Z o x f x' \ u ( x 1 ^ 
-x 2 -2- — +- — 



(16) 
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Since in the definition of fwp^^jj according to equation only G actually depends on x 
and i, the kernel G, as defined in equation (|16|) . must also fulfil the time-dependent SE. Inser- 
ting (|16p into the SE ([T|) shows that z(t) and u(t) not only fulfil the same equation of motion 
as r](t) and A(t) but, in addition, are also uniquely connected via the relation 

zu — iiz = 1. (17) 

Expressing u and z in polar coordinates according to u = a cos ip and z = a sin cp shows that this 
coupling is identical to relation (|13p that connects the amplitude and phase of A. From equa- 
tion (117p it also follows that, e.g., with the knowledge of z, the quantity u can be calculated as 

z A 



The last necessary step is to explicitly perform the integration in (|14p . using (|15|) and (|16p . to 
yield the WP solution in the form 



\V4 / 1 \ 1/2 (im \z cy {x-^z) 2 



V WP (x,t)= — ) — - exp 

\7rn/ \ u + jz / I In 



-x 



z z(u + iz) 



Comparison with the same WP, written in the form given in equation ([2]), shows that the 
relations 

777 

-i7(t) (18) 



a p 



and 



2h z 1 _ A 
m ^ z zX A 

are valid, where A = u + iz is identical to our linearization variable and equation (|17p has been 
applied. So, from equation (|18p . it then follows that the imaginary part of A is, apart from 
a constant factor, just the particle trajectory. 

The equivalence between deriving the time-dependent Green's function via a Gaussian ansatz 
or via Feynman's path integral method has been mentioned in [6], where also the relation to the 
Ermakov invariant is considered. 

In conclusion, one can say that the complex quantity A contains the particle as well as the 
wave aspects of the system. In polar coordinates, the absolute value a of A is directly connected 
with the quantum mechanical position uncertainty; in cartesian coordinates, the imaginary part 
of A is directly proportional to the classical particle trajectory r/. Absolute value and phase, or 
real and imaginary parts, of A are not independent of each other but uniquely connected via the 
conservation laws (fT3l) and (fT7|) . respectively. 



3.3 Quantum-classical connection 

Further insight into the relation between the classical and quantum dynamics of the system can 
be gained by rewriting the invariant (|10p . with the help of equation (|18p. in terms of z and z 
instead of r] and r), 



1 fa p \ 

Il = 2 l^rJ 



(za — za) 2 + ( — 
a 



const. 
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Since ^ = simp, for 7l to be constant it is necessary that (za — za) 2 = cos 2 ip = (^) 2 is valid. 
So, up to a ± sign, one obtains 

u = za 2 — zaa = a 2 ( z z ) . (19) 

V « / 

With this form of u, a certain asymmetry in the exponent of the Feynman kernel (|16p . namely | 
as coefficient of x 2 compared to as coefficient of the initial quantity x' 2 , can be explained. 

According to (fT9|) . | can be written as 

z 1 u a , . 

z or(i) z a 

This shows that, in the case when a is time-dependent, not only a 2 , must be replaced by a 2 (f), 
but also an additional term — must be taken into account. For constant a = ao or at t = 0, 
relation (1201) reduces to 



z 1 it 2 . ma . 

- = — it— or u = a n z = ry, (21) 

2 agz po 

i.e., li is simply proportional to ?). Note the explicit occurrence of ao, the initial position 
uncertainty, because it has important consequences for the time-dependence of a(t). Inserting 77, 
as given in (I2T1) , into a 2 (f ) = u 2 + z 2 yields (with v$ = ^ and (3$ = \ ) 

a 2 (t) = 4W + m=^ 2 )- (22) 

This shows that the quantum mechanical uncertainty of position (at any time t) can be expressed 
solely in terms of the classical trajectory n(t) and the corresponding velocity 77(f), if the initial 
velocity v$ and the initial position uncertainty, expressed by ao are known. 

This explains why Feynman's procedure [7] of deriving his kernel based only on the classical 
Lagrangian provides the correct time-evolution of the system since the time-dependence enters 
only via the classical variables r](t) and rj(t). However, the importance of the initial position 
uncertainty «o f° r the quantum dynamics should not be underestimated. 

The influence of the initial uncertainty becomes clear if one inserts the expressions for 77(f) 
and 77(f) into equation (|22|) . 

a) For the free motion, one obtains with 77 = 7j f, f] = vq: 

a 2 ir = ag [1 + (A)f) 2 ] ; 

b) for the HO with r\ = ^ sinwf and 7) = vq coswf, equation ([22"]) yields 

2"! 



a uo = a o { cos2 wf + I — sin ut 



Po 



u 

Only if the initial state is the ground state, is /3q = 2m (^) = = u va ^°- an d, hence, 
a = ao, i.e., the WP width remains constant; in all other cases it oscillates. This oscillating 
WP width corresponds to the general solution of the Riccati equation and yields in the limit 
u! — ► the correct spreading WP width of the free motion WP, 

Urn a| (f) = a 2 [l + (f3 t) 2 } = a 2 r (f), 
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whereas the WP usually presented, with constant width ao 
a plane wave, no WP! 

The time-derivative of a 2 , or, 



1 /2 

(^) , provides in this limit only 



■o 



d_ 

drj 



V{rj) 



respectively, only vanishes if the term in square brackets is equal to zero, which depends 
on jfrjViv)- Therefore, for V = (?) = 0), a ^ is always valid; for the HO (r/ = — u 2 rj), 
a = is only valid for (3q = \ = ui, otherwise, and in particular for uj = u(t) (which describes, 

e.g., the motion of an ion in a Paul trap 0), d / always holds. So, obviously the initial 
value ao of the position uncertainty plays an important role in the qualitative behaviour of the 
quantum aspect of the dynamics. A more detailed discussion of this problem can be found in [9], 
also including dissipative effects. 



4 Ermakov invariant and Wigner function 



In the previous section, it has been shown how the real and imaginary parts (in cartesian 
coordinates) of the complex variable A(i), that allows the linearization of the complex Riccati 
equation ([J]), enter the Feynman kernel that describes the transformation of an initial quantum 
state into a state at a later time t as time-dependent parameters. But, also in polar coordinates a 
and tp (which are not independent of each other but related via <p = ^-), X(t) is related, via the 
Ermakov invariant (|1Q|) , with another description of quantum systems that shows close similarity 
with the classical phase space description of dynamical systems, namely the Wigner function. 
In order to show this connection the invariant /l shall be written explicitly in the form 



II 



a 2 + 



Q- 



V 



2aarjf] + a 2 ?) 2 



with terms bilinear in r\ and rj. How the coefficients of these terms are related with the quantum 
uncertainties, and thus with a, cp and A, follows from 



AA* 



AA* 



a 



m 



2%i 



• 2 , 2-2 

a + a ip 



2m 2 

ir {x )l ' 



d 

— (AA*) = 2aa = 2 



— ) = t([x,PV 

yij ft 



)l = -(xp + px) L . 

h 



(23) 
(24) 
(25) 



With the help of these relations, the invariant takes the form 
d 



II 



1 

2 
1 

mh 



*J2 



XX*!] 



(XX*)r}T) + AA*?? 2 



dt 



LV 



([x,p] + ) h r)(mr)) + (x )l(?™7) 



(26) 



The connection with the Wigner function becomes obvious if one performs the Wigner trans- 
formation [10] of our WP ([2]) according to 



W(x,p,t) 



1 



2vr/i 



dq eW/"^ ( * + h *)%P (x ~ h t 
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yielding the Wigner function in the form 



■Kh { \ yi J 2h A yi h\yi 



P 2 , 2 Ate 

+ - — xp 



Using relations ([23]) - ([25]) , this can be expressed as 

1 f 2 9 

W(x,p,t) = — exp < [(P 2 )lx 2 - ([x,p] + ) L xp+ (x 2 )lP 2 ] 

In particular, at the origin of the phase space, i.e., for x = and p = 0, where x 2 
p 2 — ► (mf]) 2 and — ► 7717777, one obtains 



(27) 



?7 2 , 



1 f 2m 

^(0,0,0 = ^ exp --/ L 



const 



with II as given in ([26]) which fulfils f W = since 7l is an invariant. For 1 / and p / 0, 
W(x,p, i) takes the form (|27|) and the equation of motion is given, as expected, by 







W(x,p,t) 



p dW dV dW 
+ 



dt v y m dx dx dp ' 

i.e., a continuity equation for an incompressible medium, also called Liouville equation in phase 
space. 

Further details, also concerning the different forms to express the Ermakov invariant in 
the exponent of the Wigner function and the physical interpretation of these forms is given 
in [11] . Other attempts to construct connections between the Ermakov invariant and the Wigner 
function can be found in [6]. 



5 Riccati and Ermakov equations in time-independent 
quantum mechanics 

5.1 Factorization of the invariant and creation/annihilation operators 

Formal similarities between the Ermakov invariant and the algebraic treatment of the HO using 
creation and annihilation operators (or complex normal modes, in the classical case) and its 
generalization in the formalism of SUSY can be found if Jl is written in a form that allows for 
factorization 



1 



-a 



with 



and 



A 



A* 



a 

r\ 77 

a 



a \ .1 
-77 ] -1— 77 



1 



a 2 AA* 



a' 



a \ 1 



a 



2h 

—y } v 

777 



2H 



777 



-y v 



(28) 



(29) 



For the HO with time-independent frequency u)$ and constant WP width «o, the real part 
of (f^y) vanishes (d = 0) and the imaginary part is simply \ = ip = ujq, so equations ([28]) 

and (I29p turn into 



.4, 



= 77 - lUJQl], 



A* = i] + iu r)- 



(30) 
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These expressions are, up to constant factors, identical with the complex normal modes or, in 
the quantized form, with the creation and annihilation operators of the HO. In this context, 
it should be mentioned that the Ermakov invariant has also been quantized, where a and a 
remain c-numbers, whereas position 77 and momentum mrj are, following the rules of canonical 
quantization, replaced by the corresponding operators, i.e., rj —> q op = q, mi] — > p op = f ^ (for 
details see, e.g., [Hj). Expressions (i30j) would then turn into the operators 

1 1 

-4o,o P = —Pop - iu Q q, A) = —p p + iu q. 

For comparison, the Hamiltonian operator of the HO can be written as 

1 2 m 1 2 ( t+r- 1\ 

op ' HO = 2m~ Po P + ~2 U ° q = hw Q\ b h +2 J 

with the creation/annihilation operators 
Comparison shows that 




So, A and A* (or A op and A* p ) are generalizations where the constant factor ±iuio, in front 
of q, is replaced by the complex time-dependent functions (f^y) or (f^y*), respectively. The 
factorization of the dynamical invariant and the relation to creation/annihilation operators is 
also discussed in [6]. 



5.2 Riccati and SUSY 



A different generalization of the creation/annihilation operators is found in SUSY where, essen- 
tially, the term linear in the coordinate q is replaced by a function of q, the so-called "superpo- 
tential" W(q), leading to the operators 



1 

V2 



W(q) T i 



Pop 



In this case, the term uj$q with constant ujq is replaced by a real, position-dependent func- 
tion W{q). The operators fulfil the commutator and anti-commutator relations 



\B-,B + U 



h dW 



fm dq 

The super symmetric Hamiltonian 



{B-,B+l 



Pip 
m 



a. 



SUSY 



#i 
H 2 



can be expressed with the help of B^ in the form 

h 2 d 2 



H t = B+B~ 



2m dq 



2 +V 1 (q), 



Ho = B~B + 



^d 2 
2m dq 



2 +V 2 (q). 
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A detailed discussion of the formalism of SUSY can be found, e.g., in [13]; for our discussion, 
only the aspects mentioned in the following will be necessary. Important in this context is 
that the supersymmetric partner potentials V\{q) and V2(q) fulfil real Riccati equations, which 
follows directly from the definition of B^, i.e., 



1 

2 



V 2 



W 2 



w 2 + 



h dW 



m aq 
h dW 



in 



(32) 
(33) 



The energy spectra of Hi and H 2 are identical apart from the ground state. H\ has the 



ground state E^ = 0, whereas the ground state E^ 1 of H 2 is identical with the first excited 



(2) 



state E± of H\. The ground state wave function of Hi, ^"q 1 , has no node and determines the 
superpotential via 



(i) 



W 



n 



m q. 



(i) 
o 



From equations (132p and (133p . then, the partner potentials V\ and V 2 follow. On the other 
hand, ^q 1 ^ is connected with V\ via the solution of the equation H±^^ = 0, i.e., 



(i) 



1*0 



2m dq' 



W z 



h dw 

m dq 



(i) 



^(1)^(1) 



The connection between the spectra of H\ and H 2 , i.e. E^ and E^ \ and the corresponding 



wave functions, VEvi and 
tors according to 



is determined via the generalized creation/annihilation opera- 



(i) 

n+1 



-.B + M 2) and *1 2) 



E 



(2) 



E 



(i) 

n+1 



(1) 
n+1' 



(34) 



where B + creates a node and B annihilates a node in the wave function. So, e.g., the first 
excited state VP^ of H\ (which has one node) can be obtained from the ground state of H 2 
(which has no node) by applying B + onto it as described in (|34p . 

In order to obtain the higher eigenvalues and eigenfunctions, in the definition of W = W\ 



(2) 

must be replaced by *X , leading to W; 



(2) 



etc., i.e., 



W, 



h 



dq ^0 



(s) 




with the corresponding operators 



l 

71 



W S =F 



ft d 



m dq 



thus creating a hierarchy that provides all eigenvalues and eigenfunctions of the Hamiltonians H\ 
and H 2 . 

In the context of this paper, only two systems with exact analytic solutions shall be considered 
explicitly, namely, the one-dimensional HO (with constant frequency lo = ujq) and the Coulomb 
problem. The latter case, a three-dimensional system with spherical symmetry (V(f) = V(r) = 
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— — ), can be reduced to an essentially one-dimensional problem via separation of radial and 
angular parts. Using the ansatz <& n i m {r) = - i S n i(r)Yi m ( , &,ip) = R(r)Yi m {d , tp) for the wave 
function (with n = total quantum number, I = azimuthal quantum number, m = magnetic 
quantum number, r, $,(p = polar coordinates), the energy eigenvalues E n of the system can be 
obtained from the radial SE 

with the effective potential 

Veil = V{r) + — — = + — — . (35) 

2mr z r 2mr^ 

The superpotential W, the energy eigenvalues E n and the supersymmetric potential V\ for 
the systems under consideration are given by: 

a) HO: V(q) = ^uj 2 q 2 (eigenfunctions ^ n (q): Hermite polynomials) 

W = ujq, 

1 



E n = hw In + 2 

Vl = ^V-^ = V(q)-E . (36) 

2 

b) Coulomb potential: V(r) = — ^- (eigenfunctions ^ n i(r): Laguerre polynomials) 
y/^e 2 (/ + l)h 



W 

E„> 



(I + l)h ^Jrnr 

2 / 2\ 2 
mc I e x 



2 \hc) (n' + / + l) 2 ' 



V- e " I lil + l)h2 I ^ 1 V F f37^ 

In the second case, the radial quantum number n! occurs which indicates the number of 
nodes in the wave function and is connected with the total quantum number n, that actually 
characterizes the energy eigenvalue, via n = n' + I + 1. 

Particularly the quantities V\ , given in equations (|36|) and (|37|) , shall be compared with similar 
expressions obtained in the next subsection where a nonlinear formulation of time-independent 
quantum mechanics is presented. 



5.3 Ermakov and nonlinear time-independent quantum mechanics 

In the following discussion of a nonlinear formulation of quantum mechanics, that is essentially 
based on the work of Reinisch [2] , formal similarities with SUSY (in the time- independent case) 
and the complex Riccati formalism (in the time-dependent case) shall be pointed out. 

Starting point is Madelung's hydrodynamic formulation of quantum mechanics [H] that uses 
the polar ansatz 
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for the wave function VP(r, t), turning the linear SE (JTJ) into two coupled equations for the 
amplitude a(f) and the phase S(f,t), i.e., the Hamilton-Jacobi-type equation 

^+^ + „-|^U0 ,38) 
(with V = Nabla operator, A = Laplace operator) and the continuity equation 

4-a 2 + -V(a 2 VS) = 0, (39) 
at m 



where a 2 = ^*vj/ = g(r, t) is the usual probability density. For stationary states, the energy of the 

d , 
at' 



system is related with the action S via JtS 1 = —E = const and the density is time-independent, 



i.e., 

where it subsequently follows that the second term on the lhs of equation (|39|) must also vanish. 
In the usual textbook treatment, this is achieved by taking V5 = 0, thus turning equation (|38|) 
into the conventional time-independent linear SE 

h 2 

Aa + Va = Ea. (40) 

2m 

This is, however, not the only possibility of fulfilling equation ([39]) since, also for VS / 0, 
this can be accomplished if only 

VS = % (41) 

with constant C. In this case, equation (138j) takes the form of a nonlinear Ermakov equation 

2m, „ fC\ 2 1 



Aa + - (E - V > a = {j) ^ < 42 > 

The corresponding complex Riccati equation, equivalent to equation ([4|) in the time-depen- 
dent case, is given here by 



(43) 



where the following substitutions must be made 

Considering first the one-dimensional HO, and introducing the dimensionless variable £ via 
( = \ko\q with Hko = Po = ±V2mE, V(() = V[q(()] and a = -^a, equation pSj) acquires the 
familiar form 

" + { l ~l) a -^ (44) 

A similar formulation of the time-independent SE in terms of this equation, but within 
a different context and different applications has also been given in [TB]. In another paper [TB] 
the relation between the Ermakov equation (|44p and the time-independent SE has been extended 
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to also include magnetic field effects. The nonlinear differential equation (|44p has also been 
used to obtain numerical solutions of the time-independent SE for single and double-minimum 
potentials as well as for complex energy resonance states; for details see [T7j. Here we want to 
concentrate on the similarities between the time-dependent and time-independent situation, in 
particular with respect to SUSY, as mentioned in Sections 5.1 and 5.2. 

Following the method described in [2], from the real solution onl(C) °f this nonlinear Ermakov 
equation (|44p the complex solution aiXC) °f t ne linear SE (|40p can be obtained via 

ol(C) = «NL(C)exp j-^-sj = a N L(C)exp|-z^ dQ' ^ ^ j , (45) 

from which a real (not normalized) solution of the time-independent SE can be constructed 
according to 

5 L (C) = K[a L (C)] = \ |a N Le* s + a NL e"^l = a NL cos ( [ d£ 2 1 — ] . (46) 

2 L J Vice a NL(C')/ 

So far, the energy E occurring in equation (|44|) is still a free parameter that can take any 
value. However, solving equation (|44p numerically for arbitrary values of E leads, in general, 
to solutions onl that diverge for increasing (. Only if the energy E is appropriately tuned to 
any eigenvalue E n of equation (|40p does this divergence disappear and the integral in the cosine 
of equation ([46]) takes for £ — ► oo exactly the value |, i.e., the cosine vanishes at infinity. So, 
the quantization condition that is usually obtained from the requirement of the truncation of an 
infinite series in order to avoid divergence of the wave function is, in this case, obtained from 
the requirement of nondiverging solutions of the nonlinear Ermakov equation (|44p by variation 
of the parameter E. This has been numerically verified in the case of the one-dimensional HO 
and the Coulomb problem and there is the conjecture that this property is "universal" in the 
sense that it does not depend on the potential V (see [2]). 

For comparison with the situation in SUSY, the HO and the Coulomb problem can be written 
in the form: 

a) HO: with p = (f§) 2 , E = E n = (n+ -» /i n = {2n \ l)2 and fi( 2 = = g follows: 

d+{l-fiC 2 )a = a+ (l-^J a = d-^a = ^, (47) 



where 



U n = yw 2 ? 2 -hu(n+^\= V(q) - E n . 

b) Coulomb problem: with a(r, i?, tp) = R{r)Yi m {^,ip) the radial part can be separated and, 
with the dimensionless variable £ = \ko\r with now hko = po = ±y / 2m(—E) (E < 0), the radial 
wave function can be written as X(£) = r(£)X[r(£)], which corresponds to ^ n i( r ) i n SUSY. 
This function obeys, again, an Ermakov equation, namely 

X+[^-l)x = X + l ^X = ± (48) 



where 
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is just the effective potential V e {{ from SUSY (see equations ([35]) and (f37|) ) and E = E n = — ^ i 
with n = n' + I + 1. The coefficient of the term linear in X can again be expressed with the 
help of the potential-like expression U n i as 



e 2 *(Z + l)ft 2 mc 2 /e 2 \ 2 1 _ 
r imr^ 2 \"- c / {n +1 + 1) 

In both cases, the ground state (n = 0) wave functions are real, nodeless (n' = 0) and the 
phase does not depend on spatial variables (i.e., V5 = 0). Therefore, the rhs of equations (|4"7|) 
and (j48l) vanishes since -\ oc (V5) 2 a (similar for the Coulomb problem), i.e., the nonlinear 
Ermakov equations turn into the usual time-independent SEs. In this case, comparison shows 
that for the HO and the Coulomb problem, the potential-like terms Uq are identical with the 
corresponding V\ of SUSY. For n > and n' > 0, however, U n and U n > are different from V\ and 
describe higher excited states. In SUSY, these states can only be obtained from the hierarchy 
described in the previous subsection. Here, the price that must be paid to include these excited 
states is the nonlinearity on the rhs of equations (|47p and (|48p . 

Comparing the situation in this nonlinear formulation of time-independent quantum mechan- 
ics with SUSY and the time-dependent systems discussed in the first part of this paper, one can 
see the following similarities: 

i) Comparison with SUSY: 

The real superpotential W = ~-J= (t^) is replaced by the complex "superpotential" C(q) = 
— -^L ^^jTp^ + i^f , i.e., the ground state is replaced by the absolute value |^| of any 

eigenstate and an additional imaginary part depending on the phase ^ of the wave function 
occurs, being responsible for the non- vanishing rhs of the Ermakov equations (147j) and (|48j) . 

ii) Comparison with time- dependent SE: 

As mentioned before, the complex quantity C{q) can be compared with the time-dependent 
quantity (f^y) fulfilling the complex Riccati equation (4), C(q) = ~^=^- (^2/(0) = J' or ' 
in terms of real and imaginary parts, <-> ^ and ^ oc j^p <-> ip oc 



6 Conclusions and perspectives 

As a result of this investigation, one can state that the Ermakov invariant is a quantity of central 
importance that connects different forms for the description of the dynamics of quantum sys- 
tems, such as the time-dependent SE, the time-dependent Green's function (or Feynman kernel) 
and the time-dependent Wigner function. Unlike the classical Hamiltonian or Lagrangian, this 
invariant not only depends on the classical variables such as position and momentum, but also on 
the quantum uncertainties of these quantities contained in a and d. Therefore, the initial values 
of these quantities also play an important role for the time-evolution of the quantum system, as 
has been demonstrated in the discussion of the time-dependence of the WP width or position 
uncertainty. So, the time-evolution of a typical quantum mechanical property can be totally 
described if one only knows the classical trajectory n and the classical velocity r] (including their 
initial conditions) plus the initial position uncertainty. This traces quantum dynamics entirely 
back to the classical one plus the existence of an uncertainty principle. 

So far, the discussion of time-dependent systems included only systems where the potential 
is at most quadratic in its variables (a similar treatment of the motion in a magnetic field is 
also possible; see, e.g., [18]). This might not be as restrictive as it seems at first sight since one 
may sometimes perform canonical transformations to reduce a given Hamiltonian to a quadratic 
form [19] which has been shown explicitly by Sarlet for some polynomial Hamiltonians. To what 
extent this method can also be applied in our case requires more detailed studies. Another way 
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to extend considerably the class of Hamiltonians for which an exact invariant can be found is by 
making use of generalized canonical transformations [20J. Further generalizations of applying 
the Ermakov-type invariants also in the context of field-atom interactions, systems of coupled 
oscillators including damping and/or time-dependent masses and attempts of obtaining the 
corresponding Wigner functions can be found in [21]. For a further survey of two-dimensional 
problems in this context, see also |22j . 

Similar to the factorization of the Hamiltonian of the HO in terms of (classical) complex 
normal modes, or, (quantum mechanical) in terms of creation and annihilation operators b^, 
a factorization of the Ermakov invariant is possible that looks like a complex time-dependent 
generalization of this formalism. The major difference compared with the conventional case is 
the replacement of Hujq in front of the spatial variable q in equation (|3ip by the complex time- 
dependent quantity (f^y) , which fulfils the Riccati equation (j3|), or its conjugate complex (^y*) , 
respectively. 

Another generalization of the creation/annihilation operator formalism of the HO, concerning 
the space-dependence instead of the time-dependence, is found in SUSY where the term linear 
in the spatial coordinate q is replaced by a function of q, the so-called "superpotential" W(q). 
The generalized creation and annihilation operators are obtained by replacing the term cu^q 
in b^ by (up to constants) W(q), where W(q) fulfils the real Riccati equations (f32j) and (j33j) . 

In a nonlinear formulation of time- independent quantum mechanics proposed by Reinisch [2] , 
the amplitude of the wave function fulfils a real nonlinear (space-dependent) Ermakov equation, 
that can be connected with the time-independent SE written as a complex space-dependent 
Riccati equation (see (|43p ) together with a kind of conservation law (see (|4ip). similar to the 
conservation of "angular momentum" in the complex plane (see (I13p ) in the time-dependent 
case. In SUSY, the superpotential is initially determined by the amplitude of the ground state 
wave function (i.e., without any contribution from the phase) and the excited states can be 
obtained by some hierarchy based thereupon, in the nonlinear formulation of time-independent 
quantum mechanics, the ground state wave function is replaced by the absolute value of any 
eigenstate plus an imaginary part depending on the gradient of the phase of this state. Similar 
to the time-dependent situation, this looks like a complex generalization where, now, the real 
superpotential W(q) is replaced by a complex term that not only depends on the amplitude of 
the wave function, but, due to an additional imaginary contribution, also on (the gradient of) 
its phase. Further clarification of these facts will be subject of forthcoming studies. 
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